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1. Introduction 

In this paper, we address the peculiarities of criticality under an influence of the random 
anisotropy of structure. To be more specific, given a reference system is a 3d magnet with 
n-component order parameter which below the second order phase transition point T c 
characterizes a ferromagnetic state, what will be the impact of random anisotropy [1-3] 
on the critical dynamics [4, 5] of this transition? It appears, that contrary to the general 
believe that even weak random anisotropy destroys ferromagnetic long-range order at 
d — 3, this is true only for the isotropic random axis distribution [6]. Therefore, we will 
study a particular case, when the second order phase transition survives and, moreover, 
it remains in the random-Ising universality class [7, 8] for any n. A particular feature 
of 3d systems which belong to the random-Ising universality class is that their heat 
capacity does not diverge at T c (it is the isothermal magnetic susceptibility which 
manifests singularity) [9]. Again, general arguments state [10,11] that for such systems 
the relaxational critical dynamics of the non-conserved order parameter coupled to a 
conserved density, model C dynamics, degenerates to purely relaxation model without 
any couplings to conserved densities (model A). Nevertheless, this statement is true only 
in the asymptotics [12,13] (i.e. at T c , which in fact is never reached in experiments or in 
simulations). As we will show in the paper, common influence of two different factors: 
randomness of structure and coupling of dynamical modes leads to a rich effective critical 
behavior which possesses many new unexpected features. 

Dynamical properties of a system near the critical point are determined by the 
behavior of its slow densities. In addition to the order parameter density ip these are 
the conserved densities. Here, we consider the case of one conserved density m. For the 
description of critical dynamics the characteristic time scales for the order parameter, 
t v , and for the conserved density, t m , are introduced. Approaching the critical point, 
where the correlation length £ is infinite, they are growing accordingly to the scaling 
laws 



These power laws define the dynamical critical exponents of the order parameter, z, 
and of the conserved densities, z m . The conserved density dynamical exponents may be 
different from that of the order parameter. 

The simplest dynamical model taking into account conserved densities is model 
C, [4, 14] which contains a static coupling between non-conserved n-dimensional order 
parameter (p and scalar conserved density m. Being quite simple, the model can be 
applied to the description of different physical systems. In particular, in a lattice 
model of intermetallic alloys [15] the non-conserved order parameter corresponds to 
differences in the concentration of atoms of certain kind between the odd and even 
sublattices. It is coupled to a conserved quantity - the concentration of atoms of 
this kind in the full system. In the supercooled liquids the fraction of locally favored 
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structures is non-conserved "bond order parameter", coupled to the conserved density 
of a liquid [17]. Systems containing annealed impurities with long relaxational times [18] 
manifest certain similarity with the model C as well. 

Dynamical properties of a model with coupling to a conserved density were less 
studied numerically than those for model without any coupling to secondary densities. 
It may be the consequence of the complexity of the numerical algorithms, which turn 
out to be much slower than for the simpler model. Simulations were performed for an 
Ising antiferromagnet with conserved full magnetization and non-conserved staggered 
magnetization (i.e. the order parameter) [19] and also for an Ising magnet with conserved 
energy [20]. 

Theoretical analysis of model C critical dynamics were performed by means of 
the field-theoretical renormalization group. Critical dynamical behavior of model C in 
different regions of d — n plane was analyzed by e = 4 — d expansion in first order in 
e [14]. The results lead to speculations about the existence of an anomalous region 
for 2 < n < 4, where the order parameter is much faster than the conserved density 
and dynamic scaling is questionable. Recent two-loop calculation [21,22] corrected the 
results of Ref. [23] and showed an absence of the anomalous region 2 < n < 4. 

For the 3d model C with order parameter dimension n — 1, the conserved density 
lead to the "strong" scaling: [21,22] the dynamical exponents z and z m coincide and are 
equal 2 + a/u, where a and v are the specific heat and the correlation length critical 
exponents, correspondingly. For the Ising system {n — 1) the specific heat diverges 
and a > 0. While for a system with a < 0, that is for the physically interesting 
cases n = 2,3, the scalar density decouples from the order parameter density in the 
asymptotic region. It means that for such values of n the order parameter scales with 
the same dynamical critical exponent z as in the model A and the dynamical exponent 
of the scalar density is equal to z m = 2. The importance of the sign of a was already 
mentioned in Ref. [14]. 

A rich critical dynamical behavior has already been observed in system with 
structural disorder [18,24-27]. Interest in this case is increased by the fact that real 
materials are always characterized by some imperfection of their structure. Obviously, 
that models describing their properties should contain terms connected with structural 
disorder of certain type. For the static behavior of a system with quenched energy 
coupled disorder (e.g. dilution), the Harris criterion [28] states that disorder does not 
lead to a new static universality class if the heat capacity of the pure system does not 
diverge, that is a < 0. In appears that in diluted systems a < is always the case 
(see Ref. [9]). The conclusion about influence of coupling between order parameter and 
secondary density works also in this case. The presence of a secondary density does 
not affect the dynamical critical properties in the asymptotics [10]: order parameter 
dynamics is the same as in an appropriate model A, and z m = 2. Nevertheless, as we 
noted at the beginning, the coupling between the order parameter and the secondary 
density considerably influences the non-asymptotic critical behavior [12,13]. 

We are interested in the critical dynamics of a systems with structural disorder of 
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another type, namely, random anisotropy magnets. Their properties are described by 
the random anisotropy model (RAM) introduced in Ref. [1]. In this spin lattice model 
each spin is subjected to a local anisotropy of random orientation, which essentially is 
described by a vector and therefore is defined only for n > 1. The Hamiltonian reads: [1] 

H = - •Vr'SrSr; - D ^T^rSr) 2 , (3) 
R,R' R 

where, S = (S 1 , S n ), are n-component vectors located on the sites R of a d- 
dimensional cubic lattice, D > is an anisotropy constant, x is a random unit vector 
pointing in direction of the local anisotropy axis. The short-range interaction Jr 5 r/ is 
assumed to be ferromagnetic. 

The static critical behavior of RAM was analyzed by many theoretical and 
numerical investigations which could be compared with the critical properties of random 
anisotropy magnets found in experiments (for recent review see Ref. [3]). The results 
of this analysis bring about that random anisotropy magnets do not show a second 
order phase transition for an isotropic random axis distribution. However they possibly 
undergo a second-order phase transition for an anisotropic distribution (for references 
see reviews Refs. [2,3]). Renormalization group studies of the asymptotic [7,8,29-31] and 
non-asymptotic properties [3] of RAM corroborated such a conclusion. For example, the 
RAM with random axes distributed due to the so-called cubic distribution was shown 
within two-loop approximation to undergo a second order phase transition governed by 
the random Ising critical exponents, [3,8] as first suggested in Ref. [7]. Recently this 
result found its confirmation in a five-loop RG study [31]. The cubic distribution allows 
x to point only along one of the 2n directions of the axes ki of a (hyper) cubic lattice: [29] 

1 n 

= ^ E [ 6(n) ^ ~ **) + * (n) (* + **)], (4) 

i=l 

where 5(y) are Kronecker's deltas. 

Contrary to the static critical behavior of random anisotropy magnets their 
dynamics was less investigated. Only dynamical models for systems with isotropic 
distribution were briefly discussed in Refs. [32,33]. The critical dynamics was discussed 
within model A, Ref. [34], and the dynamical exponents were calculated. However, it 
does not give a comprehensive quantitative description since it is (i) restricted to the 
isotropic distribution of the random axis and (ii) it is performed only within the first 
non-trivial order of e = 4 — d expansion. 

The model A critical dynamics of RAM with cubic random axis distribution was 
analyzed within two-loop approximation in Ref. [35] Although the asymptotic dynamical 
properties found coincide with those of the random-site Ising model, the non-asymptotic 
behavior is strongly influenced by the presence of random anisotropy [35]. 

Beside the slow order parameter an additional slow conserved densities might be 
present, for instance the energy density. Therefore considering the non-asymptotic 
dynamical behavior of the RAM an extension to model C is of interest. Indeed, there 
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exist magnets where the distribution of the local random axes is anisotropic (e.g. the 
rare earth compounds, see Ref. [3]). 

The structure of the paper is as follows: Section [2] presents the equations defining 
the dynamical model and its Lagrangian, the renormalization is performed is Section [31 
there the asymptotic and effective dynamical critical exponents are defined. In Section 
|4] we give the expressions for the field-theoretic functions in two-loop order and the 
resulting non-asymptotic behavior is discussed. Section [5] summarizes our study. Details 
of the perturbation expansion are presented in the appendix. 



2. Model equations 

Here we consider the dynamical model for random anisotropy systems described by 
(J3]) with random axis distribution (j4j). The structure of the equations of motion for 
n-component order parameter (p and secondary density [4, 14] m is not changed by 
presence of random anisotropy 

— =-r^- + ^, , = l...n, (5) 

The order parameter relaxes and conserved density diffuses with the kinetic coefficients 
f, A correspondingly. The stochastic forces 9^, 6 m obey the Einstein relations: 

< 6 Vi 0, t)B Vj {x\ t') > = 2T5(x - x')5(t - t')5 i3 , (7) 
<9 m (x,t)6 m (x',t')> =-2XV 2 5(x-x')5(t-t')5 ij . (8) 

The disorder-dependent equilibrium effective static functional TC describing behavior of 
system in the equilibrium reads: 



H- 




iv 11 
\V<fo\ 2 +r\0 o \ 2 +^|^o| 4 -^o(*^o) 2 + ^ m o + — i^?^o I <^o 1 2 - hm } , (9) 



where D is an anisotropy constant proportional to D of Eq. (j3J), f and v depend on D 
and the coupling of the usual (f) A model. 

Integrating out the secondary density one reduces © to usual Ginzburg-Landau- 
Wilson model with random anisotropy term and new parameter v and r connected to 
the model parameters f, v, 7 and h via relations: 

r — r + jh, v = v — Sj 2 (10) 

We study the critical dynamics by applying the Bausch-Janssen- Wagner approach 
[36] of dynamical field-theoretical renormalization group (RG). In this approach, the 
critical behavior is studied on the basis of long-distance and long-time properties of 
the Lagrangian incorporating features of dynamical equations of the model. The model 
defined by expressions flSJ)-© within Bausch-Janssen- Wagner formulation [36] turns out 
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to be described by an unrenormalized Lagrangian: 



/( n n ( d 

dRdtl -f ^ 0o,i0o,*+ ^°>* ( V 
^ i=i i=i ^ 

- AV 2 ^ m + if v(p ,i<fo,i 



Am V m + rh 



<Po,i + 



E 



with auxiliary response fields (f>i{t). There are two ways to average over the disorder 
configurations for dynamics. The first way originates from statics and consists in using 
the replica trick, [37] where N replicas of the system are introduced in order to facilitate 
configurational averaging of the corresponding generating functional. Finally the limit 
iV — > has to be taken. 

However we follow the second way proposed in Ref. [33]. There it was shown that 
the replica trick is not necessary if one takes just the average of the Lagrangian with 
respect to the distribution of random variables. The Lagrangian obtained in this way is 
described by the following expression: 



C = < / dUdt (fifl 



9 +f (/i-V 2 ) J ^o-f%+^i,oj]%o%o+ 



dt 







+Am V 2 m +m [ — - AV 2 ) m + Tjm (p ifi ip ifi - 



-A7m V 2 v2i,oV ? i,o+ / dt' ^ <p ifi (t)(pi, Q (t) -^-^2<Pj,o{t')<fj,o( t ')+ 



T 2 w 
3! 



(12) 



In Eq. ( TT2|) . the bare mass is fx = f — D/n, and bare couplings are u > 0, v > 0, 
w < 0. Terms with couplings u and w are generated by averaging over configurations 
and the values of u and w are connected to the moments of distribution (jlj). Therefore 
the ratio of the two couplings has to be w/u = —n. The y-term in (1121) does not 
result from the averaging procedure but has to be included since it is generated in the 
perturbational treatment. It can be of either sign. 



3. RG functions 

We perform renormalization within minimal subtraction scheme introducing renormal- 
ization factors Z a ., ai = {{a}, {5}}, leading to the renormalized parameters {a} = 
{u, v, w, y, 7, T, A} and renormalized densities {6} = {<p, (p, m, in}. For the specific heat 
we need also an additive renormalization A, n i which leads to the function 



5^( M ,A) =//Z 2 2 /i 



dn V ^ p 



(13) 
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with the scale parameter /i and factor that renormalizes the vertex with ip 2 insertion. 
From the Z-factors one obtains the (^-functions describing the critical properties 

C(W)=-^, (14) 

Relations between the renormalization factors lead to corresponding relations between 
the ^-functions. In consequence for the description of the critical dynamics one needs 
only (^-functions of the couplings, ( Ui (ui = {u,v,w,y} for i = 1,2,3,4), the order 
parameter ( v , the auxiliary field (^-insertion ( v 2 and also function B^z. In particular, 
the ^-function of the time scale ratio 

W= T - (15) 

introduced for the description of dynamic properties is related to the above ^-functions: 
X x 

Cw = 2^ - - 7 2 ^2. (16) 

The behavior of the model parameters under renormalization is described by the 
flow equations 

^ = A... (-) 

The /3-functions for the static model parameters have the following explicit form: 

= u i (e + C v + Cu i ), (18) 

h =7(| + C^ + y^)- (19) 
The dynamic /3-function for the time scale ratio W reads 

^ ~ 2 

The asymptotic critical behavior of the system is obtained from the knowledge of 
the fixed points (FPs) of the flow equations (fT7|) . A FP {a*} = {u*, v*, w*, y*, 7*, W*} 
is defined as simultaneous zero of the /^-functions. The set of equations for the static 
fourth order couplings decouple from the other /3-functions. Thus for each of the FPs 
of the static forth order couplings {u*} one obtains two FP values of the static coupling 
between the order parameter and the conserved density 7: 

7 * 2 =0 and 7 *2 = e-2C^(K}) = « (21) 

where a and v are the heat capacity and correlation length critical exponent calculated 
at the corresponding FP {u*}. Inserting the obtained values for the static FPs into the 
/5-function (120]) one finds the corresponding FP values of the time scale ratio W. 

The stable FP accessible from the initial conditions corresponds to the critical 
point of system. A FP is stable if all eigenvalues Ui of the stability matrix df3 a Jdotj 
calculated at this FP have positive real parts. The values of u>i indicate also how fast 
the renormalized model parameters reach their fixed point values. 



Pw = W(w = W(-C v - ^ ~ 1 2 BA- (20) 
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From the structure of /3-functions we conclude, that the stability of any FP with 
respect to the parameters 7 and W is determined solely by the derivatives of the 
corresponding /^-functions: 

^ = W Uw= mV (22) 



Moreover using ([19]) we can write: 



e - 2( v 2({ Ui }) 3 2 
o; 7 = + -7 5^2 (u, A) , (23) 



2 2 
which at the FP {a*} leads to: 

^l W={Q * } = -|; for 7 * 2 = 0, (24) 

^l W=M for fVO. (25) 

Therefore, a stability with respect to parameter 7 is determined by the sign of the 
specific heat exponent a. For a system with non-diverging heat capacity (a < 0) at the 
critical point, 7* = is the stable FP. Static results report that the stable and accessible 
FP is of a random site Ising type. In this case a < 0. This leads to the conclusions that 
in the asymptotic region the secondary density decouples from the order parameter. 

The critical exponents are defined by the FP values of the ^-functions. For instance, 
the asymptotic dynamical critical exponent z is expressed at the stable FP by: 

2 = 2 + Cr(K», (26) 

with 

Cr(M) = \Cp{M) - (27) 

In similar way the dynamical critical exponent z m for the secondary density is defined 
by: 



2 + UK}, 7*), (28) 



where 



Cm{{u l },i) = -i 2 B^{{u, i }). (29) 



1 
2 

While their effective counterparts in the non asymptotic region are defined by the 
solution of flow equations (jTTI) as 

z eS = 2 + Cr(Mm,7(*),Wm (30) 

4 ff = 2 + ( m (KWb 2 M). (si) 

In the limit £ — > the effective exponents reach their asymptotic values. In the next 
section we analyze the possible scenarios of effective dynamical behavior as well as check 
the approach to the asymptotical regime. 
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4. Results 

4-1. Asymptotic properties 

The static two-loop RG functions of RAM with cubic random axis distribution in the 
minimal substraction scheme agree with the results obtained in Ref. [3] using the replica 
trick and read: 



4 2 n + 2 2 1 7 o ll(ra + 2) » 

3 U = — eu-\ — u -\ vu+-uw+yu-\ — wv — u vu 

H 3 3 3 y 3 6 18 

5 (n + 2) 2 11 2 5 2 11 2 5 



Vu u~w uW Wy y 2 u 



36 9 18 6 " 12 

-vuw wyu vuy v 2 w w 2 v , (32) 

2 6 6 9 9 ( J 



n + 8 2 2 3n+14 o lln+58 2 41 2 

= - eiH — u +2vu+-wv +yv — — v — — v u-—vu - 

6 3 12 18 18 

5 2 5 31 2 11 2 5 2 17 17 
— w v--vwy-—v w — —v y-j^V v — -vuy-—wvu, (33) 



4 2 o 2 7 3 29 2 41 2 31 2 
p w = — ewH — w +2wmH — wv+yw — w w u wu w v— 

H 3 3 y 6 9 18 18 

n+10 2 11 2 5 2 17 5n + 34 5 

— — v w — —w y~—y w — -wuy — — wvu--vwy, (34) 

36 6 12 6 18 6 



3 2 4 17 o 41 2 

Py= -zy+T^y +2yu + 2yv + 2wy+ -wv-—y - —u y- 

23 2 23 2 23 2 5n + 82 2 41 2 2 

— y u y v y w v y w y—2w v— 

6 6 6 36 y 18 

n+ 18 2 41 41 5n + 82 8 

v w vwy wyu vuy — vuw, (35) 

9 6 y 9 y 18 y 3 v; 

l 9 y 2 1 9 n + 2 ~ yu wv yv n + 2 wu wy .„„,. 

r = _ M 2 + ^ + _ w 2 + v 2 +— + — +— + fu+ — +— , 36) 

v 36 24 36 72 12 12 12 36 18 12 ' v ; 



. 1 n + 2 1 y 5 2 5 „ 2 5 2 5 2 

C (fl 2=-uH v+-w H u I2v — w — y — 

Sv 3 6 3 2 12 n + 2 6 4 y 

n + 2 5 ?/m wv yv wy 

5 vu — wu— 5 5 5 5 — . (37) 

6 3 2 2 2 2 v ; 

Here, u,v,w,y stand for the renormalized couplings. 

Given the expression for the function ^2, Eq. (I3TI) . the function /3 7 can be 
constructed via Eq. (fl9l) and the two-loop expression £^2 = n/2. 

In order to discuss the dynamical FPs it turns out to be useful to introduce the 
parameter p = W/(l + W) which maps W and its FPs into a finite region of the 
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parameter space p. Then instead of the flow equation for W the flow equation for p 
arises in (fTTj) : 

^ = &(H},7,P), (38) 

where according to (|20|) 

7, p) = p(p - l)(Cr(M, 7, P) - 7X 2 (39) 

The function £p in the above expression is obtained from Eq ff27|) using the static 
function (136]) and the two loop result for the dynamic function (calculated from 
Eq ( 1A.8I) ). We get the following two- loop expression for (r- 

u + w 2 (6 In (4/3) - 1) f(n + 2) 2 2 



Cr = _ _ + y p + — ^^-^ + ^ + r j + 

— (5u 2 + (n + 2) uv + lOtoy + 3ity + 3ww + 5if 2 + 3wy) — 
36 

+ (1-31^(4/3)) + 
Pi' ~ P ~ 3 ( " 2 + 2> In (4/3) + (1 + p) In (l - p 2 )) + u + w\ + 

7V(^)(^ + (3 + />)ln(l-rt). (40) 

The two-loop result [22] for the pure model C is recovered by setting in fT40l the couplings 
u,w,y equal to zero. While setting 7 = in (j4"0]) the result for model A with random 
anisotropy [35] is recovered. The 7 2 u, 7 2 u>, 7 2 |/-terms represent the intrinsic contribution 
of model C for random anisotropy magnets. 

There are two different ways to proceed with the numerical analysis of the 
perturbative expansions for the RG functions (132]) - (135]) . pO]) . The first one is an 
e-expansion [38] whereas the second one is the so-called fixed-dimension approach [39]. 
In the frames of the latter approach, one fixed e and solves the non-linear FP equations 
directly at the space dimension of interest (i.e. at e = 1 in our d = 3 case). Whilst 
in many problems these two ways serve as complementing ones, it appears that for 
certain cases only one of them, namely the fixed-rf approach leads to the quantitative 
description. Indeed, as it is well known by now, the e-expansion turns into the ^fe- 
expansion for the random-site Ising model and no reliable numerical estimates can be 
obtained on its basis (see [9] and references therein). As one will see below, the random- 
site Ising model behavior emerges in our problem as well, therefore we proceed within 
the fixed-d approach. 

The series for RG functions are known to diverge. Therefore to obtain reliable 
results on their basis we apply the Pade-Borel resummation procedure [40] to the static 
functions. It is performed in following way: we construct the Borel image of the resolvent 
series [41] of the initial RG function /: 

/= Yl a w (uty(vty(wt) k (yt) 1 ^ 

0<i+j+l+k<2 
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0<i+j+l+k<2 



E 



0'i,j,k,iu l v j w k y l t' l+j 
T{i+j + l) 



where / stands for one of the static RG functions fl Ui , /3 7 /7 — 7 2 n/4, a^k,! are the 
corresponding expansion coefficients given by Eqs. (!32|) - (l371) . and T(i+j + 1) is Euler's 
gamma function. Then, the Borel image is extrapolated by a rational Pade approximant 
[42] [K/L](t). Within two-loop approximation we use the diagonal approximant with 
linear denominator [1/1]. As it is known in the Pade analysis, the diagonal approximants 
ensure the best convergence of the results [42] . The resummed function is then calculated 
by an inverse Borel transform of this approximant: 



As far as the above procedure enables one to restore correct static RG flow (as sketched 
below) we do not further resum the dynamic RG function j3w- 

The analysis of the static functions fi{ Ui } at fixed dimension d = 3 brings about an 
existence of 16 FPs [3,31]. Only ten of these FPs are situated in the region of physical 
interest u > 0, v > 0, w < 0. Corresponding values of FP coordinates can be found in 



For each static FP {«*} we obtain a set of dynamical FPs with different 7* and 
p*. The FPs obtained for n = 2, 3 are listed in Tableland Table [2] correspondingly. 
Stability exponents u; 7 and lo p are given in tables as well. Here we keep the numbering 
of FPs already used in Refs. [3,8,29,31,35]. It is known from the statics that FP 
XV governs the critical behavior of RAM with cubic distribution. This FP is of the 
same origin as the FP of random-site Ising model therefore all static critical exponents 
coincide with those of the random-site Ising model. Since the specific heat exponent in 
this case is negative, the asymptotic critical dynamics is described by model A. However 
the non-asymptotic critical properties of random anisotropy magnets are different from 
the random-site Ising magnets in statics [3] as well as in dynamics [35]. Moreover, for 
the model C considered here, the non-asymptotic critical behavior differs considerably 
from that of the corresponding model A as we will se below. 

4-2. N on- asymptotic properties 

The existence of such a large number of dynamical FPs makes non-asymptotic critical 
behavior more complex as in model A. We present here results for n — 3. For n = 2 
the behavior is qualitatively similar. Solving the flow equations for different initial 
conditions we obtain different flows in the space of model parameters. The projection 
of most characteristic flows into the subspace w — y — p is presented in Fig. [TJ The open 
circles indicate genuine model C unstable FPs whereas filled circles represent model A 
unstable FPs. The filled square denotes the stable FP. 

The initial conditions for the couplings u(0),v(0),w(0),y(0) for the flows shown 
are the same as those in Refs. [3,35]. We choose 7(0) = 0.1 and p(0) = 0.6. Many 
flows are affected by the two Ising FPs Vc and X.q. Inserting the solutions of the flow 




(41) 



Ref. [3]. 
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Figure 1. Projections of flows for n = 3 in the subspace of couplings w — y — p. Open 
circles represent projections of unstable FPs with non-zero 7*. Filled circles denote 
unstable FPs with 7* = 0. The filled square shows the stable FP. See Section [4~2l for 
a more detailed description. 




Figure 2. Dependence of the order parameter effective dynamical critical exponent in 
the model C dynamics on the logarithm of flow parameter. See text for full description. 




Figure 3. Dependence of the conserved density effective dynamical critical exponent in 
the model C dynamics on the logarithm of flow parameter. See text for full description. 
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Table 1. Two- loop values for the dynamical FPs of random anisotropy magnets with 
n = 2 (model C). 
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equations into the expressions for dynamical exponents we obtain the effective exponents 
z eS and ■ The dependence of z eS on the flow parameter t corresponding to flows 1-7 
is shown in Fig. [2j Similarly Fig [3] shows this dependence for the effective exponent of 
the conserved density z^- Flow 3 is affected by both FPs and X c . Therefore the 
effective exponents demonstrate a region with values which are close to those for model 
C in the case of the Ising magnet (see curves 3 in Figs. [2] and [3]). The asymptotic values 
corresponding to the FPs Vc and Xc are indicated by the dashed line. They correspond 
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Table 2. Two-loop values for the dynamical FPs of random anisotropy magnets with 



n = 3 (model C). 
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to the values asymptotically obtained in the pure model C with n = 1, since the FPs 
Vc and Xc are of the same origin, that FP of pure model C. Curves 6 correspond to 
flows near the pure FP II. Whereas curve 7 corresponds to the flow near the cubic FP 
VIII. 

The main difference of the behavior of the effective dynamical exponent z eS in 
model C from that in model A is the appearance of curves with several peaks. The 
value of the peak appearing on the right-hand side depends on the initial condition 7(0) 
and p(0). This is demonstrated in FigfjJ 
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Figure 4. Dependence of z cS on the logarithm of flow parameter for different initial 
values 7 and p 
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Figure 5. Normalized effective dynamical critical exponents of order parameter and 
conserved density z 69 (solid line), z^f (dashed line) correspondent to flow 2. 

The effective behavior of the two dynamical critical exponents for the order 
parameter and the conserved density might be quite different as one sees comparing 
Figs. [2] and [3j However, one may ask if both exponents reach the asymptotic values 
in the same way. For this purpose we introduce a normalization of the values of the 
effective exponents by their values in the asymptotics. In particular, we introduce 
notations z** 5 = z eS /z, z^ = z^ / z m for order parameter exponent and conserved 
density exponent correspondingly. Figs. [6] show behavior of normalized exponents for 
order parameter and conserved density for flows 2 and 4 correspondingly. It illustrates 
that approach to the asymptotics for order parameter exponents and conserved density 
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exponents occurs in different way for different flows, that means for different initial 
conditions. For system with small degree of disorder (small u(0) and w(0), flow 4) the 
approach of order parameter dynamical exponent to asymptotic regime is faster than 
for the conserved density one, while for system with larger amount of disorder (flow 2) 
approach of both quantities is almost simultaneous. 

5. Conclusion 

In this paper, we have studied model C dynamics of the random anisotropy magnets 
with cubic distribution of local anisotropy axis. For this purpose two-loop dynamical 
RG function (r has been obtained. On the base of static results [3] the dependencies of 
effective critical exponents of order parameter, z eS , and conserved density, z^, on the 
flow parameter were calculated. 

The two-loop approximation adopted in our paper may be considered as certain 
compromise between what is feasible in static calculations from one side, and in dynamic 
ones form the other side. As a matter of fact, the state-of-the-art expansions of the static 
RG functions in the minimal subtraction scheme are currently available for many models 
within the five-loop accuracy [9,44] but it is not the case for the dynamic functions. 
Complexity of dynamical calculations is reflected in the current situation, when the 
results beyond two loops have been obtained for model A only. The model C even with 
no structural disorder seems to be outside present manageable problems (see the recent 
review [5]). However, there are examples which demonstrate the even in two loops 
highly accurate results for dynamical characteristics can be obtained. One of them is 
given by the critical dynamics of 4 He at the superfluid phase transition [45] . Besides, 
analysis of the two-loop static RG functions refined by resummation also brings about 
sufficiently accurate quantitative characteristics of a static critical behavior in disordered 
systems [9,46]. 

In the asymptotics the conserved density is decoupled from the order parameter 
and the dynamical critical behavior of random anisotropy model with cubic random axis 
distribution is the same as that of the random-site Ising model. Crossover occurring 
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between different FPs present in the random anisotropy model considerably influences 
the non-asymptotic critical properties. Different scenarios of dynamical critical behavior 
are observed depending of the initial values of the model parameters. The main feature 
is the presence of additional peaks on the curves for the effective dynamical critical 
exponents in comparison with the effective model A critical dynamics. 

As far as the approach to the asymptotics is very slow, the effective exponents may 
be observed in experiments and in numerical simulations. The effective exponent for the 
order parameter may take a value far away from the asymptotic one (the asymptotic 
value in in our two loop calculation is z = 2.139). The same holds for the conserved 
density effective critical exponent which may be far of its van Hove asymptotic value 
z m = 2. For example one can observe values of z eS and z^f close to those for pure Ising 
model with model C dynamics. 

This work was supported by Fonds zur Forderung der wissenschaftlichen Forschung 
under Project No. P16574 



Appendix A. Perturbation expansion 

We perform our calculations on the basis of the Lagrangian defined by f|T2|) using 
the Feynman graph technique. The propagators for this Lagrangian are shown in the 
Fig. ED 



G(k,u)8(k+k')5(cu+uj% J T fc f H{k,u)S(k+k')S(uj+Lu')8 id 

C(k,u)5(k+k')5(u+u% tj k r 'T D(k,u)5(k+k')5(u+u')5 itj 



Figure Al. Propagators for constructing Feynman graphs. G(k,ui) and H(k,uj) are 
response propagators while C(k 1 uj) and D(k,Lu) are correlation propagators. 

Response propagators G(k,ui) and H(k,u) are equal to 

G(k,u) = + f(/i + k 2 )) and H(k, u) = + \k 2 ) , (A.l) 

while the correlation propagators C(k,u) and D(k,u) are equal to 

C{k,u) =2T/\-tuj+f{li+k 2 )\ 2 and D(k, u) = 2Xk 2 /\ -tu+\k 2 \ 2 .(A.2) 

The vertices defined by Lagrangian are shown in Fig. IA2I 

We obtain an expression for the two-point vertex function T l £ by keeping the 
diagrams up to two-loop order. The result of calculations can be expressed in form: 

f k, u) = -iutl^, k, u) + f * (£, fc)f . (A.3) 

Here we introduce the correlation length = p, + jh, u, v, w, y), which is defined by 

(A.4) 



dk 2 



k 2 =0 
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TA5(k + k' + k" + k"')5{u + u/ + uj" + uj"') f 7<J(fc + k' + k")S(uj+uj'+uj")5 i;j 




T 2 B5(k + k')5(k" + k"')5(uj + u/)S(u" + uj'") \X^k 2 5{k + k 1 + k")5{uj+uj'+uj")8 i j 



Figure A2. Vertices for our model. In vertex a, A stands for 
wo/3! (5i t j5i lW , + Si t i5j tm + Si, m 5j,i)/3 or yo/3\ Sijdj^S^m. In vertex b, B stands for 
Uq/3! Sij5i tm or Wq/3\ 5i jSj i8i !m . Vertices c and d originate from the coupling to the 
conserved density. 

The function T w is the static two-loop vertex function of the disordered magnet. The 
structure (1A.3I) of the dynamic vertex function of pure model C was obtained in Ref. [43] 
up to two-loop order. 

We can express two-loop dynamical function Cl$<p in the following form: 

k,u) = l+ ^(£, k, uj) + ^(£, k, uj), (A.5) 

where the one loop contribution reads: 

u + f 1 

Jk> ( — iuj + f (£- 2 + k' 2 ))(£- 2 +k> 2 ) 
while the two-loop contribution is of the form: 

n^ce, k u) = t^v 2 + 1 + k, uj) 



f (^p* + m 2 Cg 3) (L k, UJ) + f^S^, k, UJ) + 



T(—vu + — + — + —)W^ k, uj)+T(- + y + 2—) x 



f(^)f (wg^ce, ^) + ^) + 2wg^ (fj ^ w) ) , (A . 7) 

In (1A.6I) and (1A.7[) the expressions for the integrals Jc, W^^, C*( T3 ) and 5 of the 
pure model C are given in the Appendix A.l in Ref. [22], while the contributions for 
\ty(CDt) are p resen ted in the Appendix of Ref. [13]. 

Following the renormalization procedure for we obtain the two-loop 

renormalizating factor Z^: 



_,2u + w 7 2 W 1 
Z *- 1+ ~e^~eTTw + 7 2 



1 1 + W 1 l + W V2 
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4 2 9 n + 2 yu vw 



fn + 2 \\ 2 W / W \u + w 2 W 2 2 




1 / , x 2 W \ 1 /n+2 2 ^ |, 2 \ / 4 1\ 

-vw— (u + w)7 > — v H m — 

4 v M 1+W J £ \ 12 6 4 y V 3 12 / 

Y — In ln(l + W) \nW 

e 3 ' 1 + W 2 1 + W 2 v ; 21 + W 



(A.8) 
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